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Abstract 

The  innovations  or  observation  minus  forecast  (O-F)  residuals  produced  by  a  data 
assimilation  system  provide  a  convenient  metric  of  evaluating  global  analyses.  In  this 
study,  O-F  statistics  from  the  Global  Ozone  Assimilation  Testing  System  (GOATS)  are 
used  to  examine  how  ozone  assimilation  products  and  their  associated  O-F  statistics 
depend  on  input  data  biases  and  ozone  photochemistry  parameterizations  (OPP).  All 
the  GOATS  results  shown  are  based  on  a  6-h  forecast  and  analysis  cycle  using  obser¬ 
vations  from  SBUV/2  (Solar  Backscatter  Ultraviolet  instrument-2)  during  September- 
October  2002.  Results  show  that  zonal  mean  ozone  analyses  are  more  independent 
of  observation  biases  and  drifts  when  using  an  OPP,  while  the  mean  ozone  O-Fs  are 
more  sensitive  to  observation  drifts  when  using  an  OPP.  In  addition,  SD  O-Fs  (stan¬ 
dard  deviations)  are  reduced  in  the  upper  stratosphere  when  using  an  OPP  due  to  a 
reduction  of  forecast  model  noise  and  to  increased  covariance  between  the  forecast 
model  and  the  observations.  Experiments  that  changed  the  OPP  reference  state  to 
match  the  observations  by  using  an  “adaptive”  OPP  scheme  reduced  the  mean  ozone 
O-Fs  at  the  expense  of  zonal  mean  ozone  analyses  being  more  susceptible  to  data 
biases  and  drifts.  Additional  experiments  showed  that  the  upper  boundary  of  the  ozone 
DAS  can  affect  the  quality  of  the  ozone  analysis  and  therefore  should  be  placed  well 
above  (at  least  a  scale  height)  the  region  of  interest. 


1  Introduction 

Ozone  is  a  radiatively  important  trace  gas  in  the  middle  atmosphere.  As  global  nu¬ 
merical  weather  prediction  (NWP)  models  extend  upward  and  become  more  reliant  on 
satellite  observations,  accurate  modeling  of  the  radiative  environment,  including  ozone, 
seen  by  these  models  becomes  more  important  for  optimal  use  of  satellite  radiances 
(e.g.,  Derber  and  Wu,  1998;  John  and  Buehler,  2004). 

A  global  ozone  data  assimilation  system  combines  a  global  forecast  model  with 
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ozone  observations  to  produce  a  global  gridded  ozone  analysis.  The  ozone  forecast 
can  fill  in  where  observations  are  lacking,  and  aids  in  data  quality  control  by  provid¬ 
ing  a  realistic  background  field.  The  observations  can  correct  forecast  errors,  thereby 
providing  accurate  ozone  initial  conditions  for  the  next  forecast  and  analysis  cycle.  In 
addition  to  improving  derived  forecast  products  such  as  surface  ultraviolet  (UV)  indices 
(  ong  et  al ,  1996),  ozone  assimilation  and  forecasting  improve  model  stratospheric 
heating  rates  (and  hence  temperatures)  and  can  be  used  to  diagnose  transport  in  the 
lower  stratosphere,  thereby  improving  wind  assimilation  (  ,  00;  Jane 

2003).  These  and  other  motivations  for  ozone  assimilation,  as  well  as  previous  work  on 
the  problem,  are  summarized  by  Rood  (200f ).  More  recent  work  has  focused  on  the 
assimilation  of  new  stratospheric  ozone  observations  to  improve  the  ozone  analysis 
(e.g.,  Stajner  and  Wargan,  2004;  Wargan  et  al.,  2005;  Stajner  et  al.,  2006),  and  gen¬ 
eration  of  multidecadal  ozone  reanalysis  fields  (Dethof  and  Holrr ,  2004;  Oikon' 
and  O’Neill,  2006).  Ozone  data  assimilation  has  also  been  used  recently  to  evaluate 
different  parameterizations  of  stratospheric  ozone  photochemistry  (Geer  et  al.,  2006b). 

In  improving  stratospheric  assimilations  further,  Rood  (2005)  concludes  that  details 
of  the  various  assimilation  algorithms  are  now  secondary  considerations  compared 
to  the  more  general  issues  of  identifying  and  eliminating  model/data  biases  and  (re- 
latedly)  improving  physical  parameterizations  (see  also  Povalarapu  et  al.,  2005).  We 
focus  on  those  broader  issues  here  in  a  series  of  experiments  with  an  ozone  assim¬ 
ilation  system.  Specifically,  we  investigate  how  the  assimilation  system  responds  to 
improvements  in  ozone  “physics”  by  comparing  results  between  experiments  with  and 
without  a  parameterization  of  ozone  photochemistry  in  the  forecast  model.  Output  from 
each  assimilation  run  is  analyzed  objectively  using  standard  output  parameters,  such 
as  mean  and  standard  deviation  (SD)  of  the  innovations  (observation  minus  forecast 
residuals).  Additional  experiments  with  and  without  parameterized  ozone  photochem¬ 
istry  examine  how  the  mean  innovations  (i.e.,  the  bias  between  the  ozone  forecast 
model  and  ozone  observations)  respond  when  fed  with  unbiased,  then  biased  versions 
of  the  same  set  of  satellite  ozone  observations. 
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We  begin  by  describing  our  ozone  assimilation  system  (Sect.  2),  including  a  brief 
summary  of  the  observations  (Sect.  ),  and  assimilation  experiments  (Sect.  ).  We 
present  results  of  the  experiments  in  Sect.  3  and  the  discussion  of  the  results  in  Sect, 
followed  by  the  conclusions  in  Sect.  . 


2  Global  Ozone  Assimilation  Testing  System 

To  allow  us  to  develop  and  refine  an  ozone  assimilation  and  forecasting  capability, 
we  have  developed  a  Global  Ozone  Assimilation  Testing  System  (GOATS),  depicted 
schematically  in  Fig.  .  GOATS  couples  the  Goddard  Earth  Observing  System  (GEOS) 
ozone  data  assimilation  system  (DAS)  with  the  NOGAPS-ALPHA  (Navy  Operational 
Global  Atmospheric  Prediction  System  -  Advanced  Level  Physics  with  High  Altitude) 
global  forecast  model.  The  components  are  discussed  in  more  depth  in  Sects.  and 
2.2,  respectively. 

To  start  each  6-h  analysis  cycle,  the  NOGAPS-ALPHA  General  Circulation  Model 
(GCM)  is  initialized  with  archived  meteorological  fields  produced  from  the  operational 
NOGAPS  using  the  Naval  Research  Laboratory  (NRL)  Atmospheric  Variational  Data 
Assimilation  System  (NAVDAS).  This  approach  simplifies  the  ozone  assimilation  by 
avoiding  a  complete  meteorological  assimilation  for  each  model  run.  This  simplification 
offers  two  advantages.  First,  it  decouples  the  meteorological  assimilation  from  the 
ozone  assimilation,  enabling  us  to  focus  on  the  parameters  relevant  to  ozone.  Second, 
it  more  easily  facilitates  the  integration  of  the  GEOS  ozone  DAS  into  the  GOATS. 

GOATS  is  similar  to  the  off-line  ozone  data  assimilation  systems  (see,  e.g., 
et  al.,  2001)  in  that  the  meteorological  analysis  is  independent  of  the  ozone  assimi¬ 
lation,  however,  GOATS  follows  the  approach  of  Stajner  et  ai.  (2006),  where  the  full 
GCM,  rather  than  an  ozone  transport  model,  is  used  to  dynamically  advect  ozone  to 
the  next  analysis  time.  This  means  that  GOATS  has  some  advantages  as  well  as  disad¬ 
vantages  over  a  typical  off-line  ozone  DAS.  Advantages  include  the  ability  to  influence 
model  dynamics  by  changing  shortwave  heating  and  longwave  cooling  due  to  ozone. 
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In  all  the  GOATS  results  shown  here,  this  interaction  of  ozone  with  dynamics  through 
shortwave  radiative  heating  and  longwave  cooling  was  activated,  however  the  effects 
are  not  discussed  in  this  paper.  Another  advantage  for  GOATS  is  that  time  interpola¬ 
tion  of  the  meteorological  analyses  is  not  necessary,  as  the  dynamical  model  predicts 
the  meteorological  fields  along  with  the  ozone  at  each  model  time  step.  The  potential 
improvements  in  constituent  transport  from  using  GCM  winds  are  discussed  by  Roc 
(2005). 

A  disadvantage  for  GOATS  is  that  the  dynamical  model  must  be  initialized  at  the 
start  of  each  analysis  cycle.  This  is  an  important  consideration  for  GOATS  because 
the  current  NOGAPS  operational  analyses  are  capped  at  10  hPa.  As  shown  schemati¬ 
cally  in  Fig.  ,  we  circumvent  this  restriction  in  the  GOATS  by  using  GEOS4  (Goddard 
Earth  Observing  System  version  4)  meteorological  analyses  from  1 0-0.2  hPa  or  from 
10-1.0hPa.  Above  the  top  analysis  level,  meteorological  fields  from  the  top  analysis 
level  are  blended  with  climatological  fields  from  (JARS  (Upper  Atmosphere  Research 
Satellite;  Swinbank  and  Ortland,  2003)  and  CIRA  (COSPAR  International  Reference 
Atmosphere;  Fleming  et  al.,  1990)  as  described  in  Eckermann  et  al.  (2004).  These 
global  meteorological  fields  are  then  balanced  using  the  NOGAPS-ALPFIA  Nonlinear 
Normal  Mode  Initialization  (NNMI)  scheme,  before  starting  the  forecast.  Section 
examines  the  effects  of  lowering  the  top  GEOS4  analysis  level  from  0.2  hPa  to  1.0  hPa 
and  eliminating  the  NNMI. 

Figure  illustrates  a  GOATS  6-h  analysis  cycle.  The  initial  meteorological  fields  are 
set  as  described  above,  while  ozone  is  taken  from  the  previous  GOATS  ozone  analysis 
(or  from  zonal-mean  climatology  for  the  first  cycle).  The  fields  are  advanced  6  h  using 
the  NOGAPS-ALPFIA  GCM.  The  6-h  ozone  forecast  is  saved  on  the  model’s  grid  at  all 
levels  using  the  spectral  coefficients.  An  offline  code  transforms  the  data  from  spectral 
to  real  space,  providing  ozone  mixing  ratios  on  a  Gaussian  latitude-longitude  grid  as  the 
“background”  fields  for  the  data  assimilation.  These  background  fields  are  then  read 
into  the  GEOS  ozone  DAS,  where  they  are  statistically  combined  with  all  SBUV/2  (Solar 
Backscatter  Ultraviolet  Radiometer-2)  data  that  occur  within  a  6-h  window  centered  on 
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the  analysis  time:  these  observations  are  discussed  further  in  Sect.  .  The  resulting 
analyzed  ozone  fields  are  transformed  back  into  spectral  coefficients  that  provide  the 
initial  ozone  field  for  the  next  6-h  cycle,  while  the  meteorological  fields  for  the  next  cycle 
are  again  provided  from  archived  NAVDAS  and  GEOS4  analyses.  Note  that  archived 
GEOS  ozone  analyses  are  not  used  in  GOATS. 

Since  we  start  GOATS  with  zonal-mean  climatological  ozone,  it  is  necessary  to  de¬ 
termine  the  spin-up  time  required  for  this  dependence  on  the  zonal-mean  ozone  initial 
condition  to  be  eliminated.  Several  tests  compared  assimilation  runs  initialized  with 
both  climatology  and  with  a  3-D  analysis  provided  by  the  operational  GEOS  ozone 
DAS.  As  in  Geer  et  al.  (2006a),  we  found  that  ten  to  fifteen  days  of  spin-up  time  were 
sufficient  to  eliminate  most  of  the  dependence  on  initial  conditions.  When  showing  time 
average  fields,  we  begin  averaging  1 0  days  after  the  GOATS  initial  time  to  suppress  any 
significant  influence  from  the  initial  condition. 

Output  from  the  GOATS  runs  includes  the  ozone  analysis,  6-hourly  ozone  forecast, 
and  the  standard  diagnostics  provided  by  the  GEOS  ozone  DAS.  The  latter  include 
differences  between  the  incoming  observations  and  the  background  forecast  (O-F,  or 
innovations).  The  mean  and  standard  deviation  of  the  O-F  are  used  extensively  to 
monitor  the  agreement  between  model  and  observations  (  al.,  2  ). 

2.1  GEOS  ozone  DAS 

The  statistical  analysis  system,  which  is  used  to  assimilate  ozone  observations  within 
the  GOATS,  was  developed  at  NASA’s  Goddard  Space  Flight  Center  as  part  of  the 
Goddard  Earth  Observing  System  Data  Assimilation  System  (GEOS-DAS).  This  three- 
dimensional  (3-D)  ozone  assimilation  system  has  been  used  in  a  wide  variety  of 
applications,  assimilating  data  from  nadir-sounding  instruments:  SBUV/2  and  Total 
Ozone  Mapping  Spectrometer  (TOMS);  and  limb-sounding  instruments:  Microwave 
Limb  Sounder  (MLS),  Michelson  Interferometer  for  Passive  Atmospheric  Sounding  (MI- 
PAS),  Polar  Ozone  and  Aerosol  Monitor  (POAM),  and  Improved  Limb  Atmospheric 
Spectrometer-1 1  (ILAS-II).  Detailed  descriptions  of  the  assimilation  system,  error  mod- 
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eling,  and  applications  are  provided  in  Riishogaard  et  al.  (2000),  Stajner  et  al.  (2001, 
2004),  Stajner  and  Wargan  (2004),  Wargan  et  al.  (2005),  and  Stajner  et  al.  (2006).  The 

GEOS  ozone  DAS  combines  ozone  observations  with  a  background  3-D  ozone  field 
from  a  forecast  model,  using  the  Physical-space  Statistical  Analysis  Scheme  (PSAS) 
described  by  Cohn  et  al.  (1998). 

2.2  NOGAPS-ALPHA 

Operational  global  NWP  at  the  Fleet  Numerical  Meteorology  and  Oceanography  Center 
(FNMOC)  is  provided  by  the  Naval  Research  Laboratory’s  (NRL)  Navy  Operational 
Global  Atmospheric  Prediction  System  (NOGAPS:  Hogan  and  Rosmonci,  199  ).  The 
NOGAPS  Eulerian  spectral  model  currently  uses  a  T239L30  formulation  operationally, 
and  these  30  o  levels  are  depicted  in  Fig.  ,  both  the  interface  (half)  levels  (panel  a) 
and  the  profile  of  corresponding  layer  thicknesses  (black  curve  in  panel  c). 

A  comparison  of  Figs,  a  and  d  shows  that  ozone  mixing  ratios  peak  at  the  up¬ 
permost  thick  diffused  layers  of  this  L30  model,  which  makes  it  unsuitable  for  ozone 
assimilation.  Thus,  GOATS  uses  an  advanced-level  physics,  high-altitude  (ALPHA) 
nonoperational  prototype  of  the  NOGAPS  spectral  model,  known  as  NOGAPS-ALPHA. 
Only  those  aspects  of  NOGAPS-ALPHA  salient  to  GOATS  are  described  in  what  fol¬ 
lows,  since  more  detailed  descriptions  of  NOGAPS-ALPHA  have  been  provided  in  a 
number  of  recent  publications  (Eckermann  et  al.,  2004,  2006;  McCormack  et  al.,  2004; 
Allen  et  al.,  2C  ). 

2.2.1  Model  resolution 

We  typically  run  NOGAPS-ALPHA  using  an  L60  hybrid  o-p  formulation  whose  layer 
thicknesses  are  shown  in  Fig.  c  (see  also  ,  ).  However,  the 

GEOS  ozone  DAS  was  developed  and  tuned  within  the  framework  of  the  specific  L55 
hybrid  o-p  levels  of  the  Version  4  Goddard  Earth  Observing  System  (GEOS4)  model 
(Bloom  et  al.,  2005).  Thus,  to  simplify  interfacing  the  GEOS  ozone  DAS  to  NOGAPS- 
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ALPHA,  we  built  these  same  L55  o-p  levels  into  NOGAPS-ALPHA,  with  some  minor 
top  level  differences  due  to  NOGAPS  use  of  a  finite  top  level  pressure  pTOP  within  its 
hybrid  coordinate  formulation  (see  Appendix  A).  Those  L55  NOGAPS-ALPHA  half  lev¬ 
els  are  shown  in  Fig.  b,  and  the  corresponding  layer  thicknesses  in  Fig.  c.  These  L55 
levels  span  the  full  depth  of  the  stratospheric  ozone  layer  with  good  vertical  resolution 
throughout. 

To  reduce  the  computational  burden,  in  this  test  environment  we  ran  the  NOGAPS- 
ALPHA  at  a  T79  spectral  resolution  in  all  the  GOATS  experiments  reported  here. 

2.2.2  Prognostic  ozone 


NOGAPS-ALPHA  initializes  a  three-dimensional  global  ozone  mixing  ratio  field  r  that  is 
subsequently  advected  spectrally  and  can  also  be  photochemically  updated  (see,  e.g., 
McCormack  et  a  ,  2004,  2006).  In  the  GOATS,  the  prognostic  ozone  is  initialized  with 
the  latest  assimilated  ozone  at  the  start  of  each  6  h  forecast  (see  Fig.  ). 

In  common  with  other  NWP  ozone  DASs,  for  photochemistry  NOGAPS-ALPHA  uses 
a  linearized  parameterization  of  rates  of  gas-phase  ozone  photochemical  production  P 
and  loss  L,  formulated  mathematically  (Cariolle  and  Deque,  1986)  as 


(P-L)  =  (P-L)0  + 


d(P-L) 

dr 


d(P-L) 

dT 

d(P-L) 

dl 


(r  -  ra) 


(T  ~  T0) 


(2  -  10)  ■ 


(1) 


Here  r,  T  and  I  are  the  current  model  values  of  ozone  mixing  ratio,  temperature  and 
ozone  column  density,  respectively.  The  “o”  subscripts  denote  values  at  a  reference 
(or  equilibrium)  state.  The  reference  states  r0,  T0  and  I0  are  specified  by  zonal-mean 
climatologies,  while  the  four  reference  state  photochemical  coefficients  in  (1)  are  pro¬ 
vided  as  zonal-mean  lookup  tables  by  the  photochemistry  scheme. 
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NOGAPS-ALPHA  incorporates  and  can  use  any  one  of  the  various  linearized  ozone 
photochemistry  schemes  currently  in  use  by  the  world’s  major  NWP/DA  centers  and 
some  climate  models.  Preliminary  assessments  and  intercomparisons  of  some  of 
these  schemes  in  NWP/DA  applications  have  been  provided  by  McCormack  et  al. 
(2004),  Geer  et  al.  (2006b)  and  McCormack  et  al.  (2006).  A  major  motivation  for  our 
development  of  GOATS  is  to  objectively  benchmark  the  performance  of  photochemistry 
schemes  in  a  realistic  high-altitude  prototype  ozone  DA  system,  prior  to  implementa¬ 
tion  of  ozone  in  operational  DA  systems,  such  as  the  NRL  Atmospheric  Variational 
Data  Assimilation  System  (NAVDAS:  ,  ). 

NRL  has  developed  a  scheme  based  on  output  from  their  zonally-averaged  CHEM2D 
model  (  cCormack  and  Siskind,  2002).  This  CHEM2D-OPP  (Ozone  Photochemistry 
Parameterization)  scheme  is  described  by  McCormack  et  al.  (2006),  and  has  been  im¬ 
plemented  operationally  in  the  National  Centers  for  Environmental  Prediction  Global 
Forecast  System  (NCEP  GFS)  and  transitioned  to  FNMOC  for  the  operational  NO¬ 
GAPS.  CHEM2D-OPP  was  used  for  all  the  GOATS  results  with  parameterized  ozone 
photochemistry  presented  in  Sect.  . 

In  the  present  study,  the  GOATS  assimilates  SBUV/2  observations  during  the  pe¬ 
riod  of  September-October  2002.  An  examination  of  the  individual  terms  in  the  pho¬ 
tochemistry  parameterization  shows  that  the  second  term  on  the  right  hand  side  of 
Eq.  (  )  dominates  the  photochemistry  component  of  NOGAPS-ALPHA  at  altitudes 
above  lOhPa  during  this  period.  This  term  acts  as  a  photochemical  relaxation  rate, 

t0s=-  [  a(Pg~L)  |  ]  (see  Cariolle  am  Deque,  1  36;  McCormack  et  al.,  2  ),  that 
draws  the  forecasted  ozone  mixing  ratio,  r,  toward  the  model’s  specified  ozone  refer¬ 
ence  state,  rQ.  The  CHEM2D-OPP  photochemical  relaxation  rates  for  September  are 
shown  in  Fig.  b. 

In  NOGAPS-ALPHA,  the  ozone  reference  state,  r0 ,  is  specified  using  an  amalgam  of 
the  (  )  zonal-mean  climatology  at  altitudes  below  0.3  hPa,  and 

long-term  diurnally-averaged  output  from  the  CHEM2D  model  above  (see  McCormack 
et  al.,  2006,  for  details).  Note  that  the  Fortuin-Kelder  climatology  is  heavily  weighted 
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with  SBUV  data,  and  thus  is  a  good  match  for  our  GOATS  SBUV/2  assimilation  ex¬ 
periments.  This  NOGAPS-ALPHA  zonal-mean  ozone  reference  state  for  September  is 
shown  in  Fig.  a. 

2.3  SBUV/2  Observations 

All  the  GOATS  experiments  that  follow  assimilate  SBUV/2  ozone  observations.  The 
SBUV/2  instrument  and  the  inverse  model  used  to  retrieve  layer  ozone  values  are 
described  in  Bhartia  et  al  (1 996).  The  GEOS  ozone  DAS  was  developed  using  SBUV/2 
observations  (  ishogaard  et  al.,  2000;  Stajner  et  ,  2001),  and  use  of  the  innovation 
statistics  from  SBUV/2  ozone  assimilation  to  monitor  the  DAS  is  described  in  Stajner 
et  al.  (  ).  Following  (  ),  we  have  used  only  SBUV/2  layers  3-12 

in  the  GOATS  along  with  the  SBUV/2  total  ozone  column  data.  The  pressure  intervals 
of  the  SBUV/2  layers  3-12  are  given  in  Table  (see  also  Fig.  3).  These  SBUV/2  layers 
are  similar  to  the  Umkehr  layers  used  for  ground-based  ozone  observations. 

As  detailed  further  in  Sect.  1,  we  used  two  operationally  issued  versions  of 
SBUV/2  observations  with  different  biases  in  the  upper  stratosphere.  These  were  ver¬ 
sions  n16v61608  and  n16v61814.  The  GOATS  was  run  using  each  version  in  separate 
experiments  to  examine  the  effect  of  the  bias  on  the  ozone  assimilation. 

2.4  GOATS  experiments 

All  the  GOATS  experiments  shown  here  started  on  1  September  2002,  00:00  UTC  and 
ran  with  a  6-h  ozone  analysis  cycle  for  up  to  2  months.  The  global  gridded  ozone 
forecast  and  analysis  fields  were  archived  every  6h  along  with  the  O-F  residuals.  A 
reduction  in  the  O-Fs  is  a  standard  measure  of  improvement  in  a  data  assimilation 
system,  and  therefore  the  O-Fs  will  be  examined  in  detail  in  Sect.  . 

The  O-Fs  can  be  examined  globally  or  within  a  more  restricted  domain.  For  each 
level  and  region  (usually  global)  we  have  computed  an  average  (mean)  O-F,  or  bias, 
and  after  subtracting  the  average  from  each  value  in  the  region,  we  have  computed  a 


1076 


standard  deviation  (SD)  O-F  over  the  same  region. 

The  time  period  that  we  are  examining  was  highlighted  by  the  first  observed  major 
stratospheric  warming  in  the  Southern  Hemisphere  (Newman  and  Nash,  2005;  Allen 
et  al. ,  ).  Ozone  assimilation  products  during  this  time  period  have  been  examined 

by  Geer  et  al  (2006a).  Figure  shows  total  ozone  comparisons  between  the  GOATS 
T79L55  ozone  analyses  (with  CHEM2D-OPP  activated)  and  TOMS  (Total  Ozone  Map¬ 
ping  Spectrometer)  during  the  warming  event.  The  GOATS  analyses  capture  the  total 
ozone  morphology  of  this  period  including  the  splitting  of  the  ozone  hole  into  two  pieces 
on  the  24th,  the  return  to  a  single  smaller  ozone  hole  on  the  30th,  and  the  ridges  of 
high  ozone  surrounding  the  ozone  hole.  These  maps  agree  well  with  the  SBUV/2  to¬ 
tal  ozone  maps  analyzed  for  this  period  by  Kondragunta  et  al.  (2005).  Note  that  the 
GOATS  low  total  ozone  values  (~150  DU)  on  20  September  2002  are  produced  by  the 
assimilated  SBUV/2  observations  that  often  see  into  the  disturbed  low  ozone  region 
during  the  first  20  days  of  the  assimilation. 

In  the  next  section  we  provide  more  quantitative  assessments  of  GOATS  perfor¬ 
mance.  Examining  GOATS  analyses  during  this  highly  dynamic  period  will  provide  a 
stringent  test  of  the  quality  of  the  analyses. 


3  Results 

This  section  describes  some  of  the  results  from  experiments  performed  with  the 
GOATS.  We  first  examine  how  the  mean  ozone  O-Fs  depend  on  parameterized  chem¬ 
istry  and  observation  biases.  Then  we  examine  how  the  SD  O-Fs  can  be  reduced 
with  the  addition  of  parameterized  chemistry.  Next  we  examine  how  biases  can  affect 
photochemical  parameterizations  in  which  reference  climatologies  are  updated  adap¬ 
tively  using  the  assimilated  fields  from  earlier  times.  Finally,  we  examine  how  improved 
mesospheric  dynamics  improve  the  ozone  assimilation. 


1077 


3.1  O-Fs  and  parameterized  chemistry 

As  shown  by  (2  ),  the  observation-minus-forecast  (O-F)  innovations 

can  be  useful  for  diagnosing  the  performance  of  an  operational  ozone  DAS.  Here, 
we  examine  the  response  of  the  O-Fs  to  changes  in  the  forecast  model,  beginning 
with  changes  occurring  when  parameterized  ozone  photochemistry  (CHEM2D-OPP) 
is  added  to  the  forecast  model. 

Figure  shows  global  average  O-F  times  series  for  SBUV/2  layer  8  (1 .98-3.96  hPa) 
for  two  experiments:  one  with  and  one  without  parameterized  ozone  photochemistry. 
Since  t03~1  day  in  layer  8  (Fig.  3b),  one  would  expect  model  chemistry  to  play  a  role 
during  each  6-h  GOATS  cycle.  As  stated  in  Sect.  4,  at  each  analysis  time  the  6-h 
mean  of  the  O-Fs  in  Fig.  b  has  been  subtracted  before  calculating  the  standard  devia¬ 
tion  values  shown  in  Fig.  a.  Separating  the  O-Fs  into  a  mean  and  a  standard  deviation 
is  useful  because,  as  will  be  shown,  each  quantity  responds  to  model  changes  in  a  dif¬ 
ferent  way.  While  parameterized  chemistry  has  reduced  the  SBUV/2  layer  8  global 
standard  deviation  of  the  O-Fs  at  all  times,  the  mean  of  the  O-Fs  shows  slightly 
larger  absolute  values  for  the  parameterized  chemistry  case  and  the  mean  O-Fs  are 
increasing  towards  the  end  of  the  integration,  which  indicates  increasing  disagreement 
between  the  observations  and  the  model  forecast.  Compared  with  a  typical  value  for 
the  layer  of  ~10  DU  (see  Table  ),  the  differences  between  the  mean  O-Fs  in  the  two 
cases  are  relatively  small. 

3.1.1  Mean  O-F 

This  section  examines  how  parameterized  chemistry  and  data  biases  affect  the  mean 
O-F.  We  have  found  that  early  versions  of  the  SBUV/2  data  were  biased  roughly  10% 
high  in  the  upper  SBUV/2  layers.  For  example,  Fig.  shows  layer  8  ozone  as  a  function 
of  latitude  from  two  versions  of  SBUV/2  for  1 5  September  2002,  along  with  zonal-mean 
September  values  from  the  rortuin  and  Kelder  (  998)  climatology,  which  NOGAPS- 
ALPHA  uses  for  r0  in  Eq.  (  ).  The  latitudinal  structures  in  the  two  SBUV/2  data  sets  are 
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very  similar  except  for  an  offset,  with  n16v61814  matching  closely  with  the  climatology 
and  n16v61608  biased  high.  We  ran  the  GOATS  in  separate  experiments  using  each 
data  set.  The  high  bias  of  the  earlier  version  allows  us  to  test  how  data  biases  affect 
the  mean  O-Fs. 

Figure  shows  the  mean  layer  8  O-F  time  series  for  GOATS  runs  using  the  different 
SBUV/2  data  sets.  The  O-Fs  with  chemistry  (Fig.  a)  are  large  for  the  n16v61814  ob¬ 
servations  that  are  biased  high  with  respect  to  the  model  ozone  climatology  in  Fig.  . 
However,  without  chemistry  (Fig.  b)  the  mean  O-Fs  converge  to  near  zero,  indepen¬ 
dent  of  the  input  data  bias.  These  latter  results  make  sense  for  global  averaged  fields 
because  the  forecast  model  without  chemistry  can  only  change  ozone  mixing  ratios  by 
advection.  In  a  global  average  there  is  no  net  change  due  to  horizontal  advection  and 
the  vertical  mass  flux  must  average  to  zero.  It  is  still  possible  to  change  global  average 
ozone  if  the  vertical  ozone  gradients  are  correlated  with  the  vertical  motion  (a  non¬ 
zero  global-mean  vertical  eddy  ozone  flux  gradient),  however,  these  changes  are  often 
small,  as  will  be  seen  in  the  vertical  profile  plots  below.  Thus,  in  experiments  without 
ozone  chemistry,  the  forecast  model  for  global  averaged  ozone  in  a  given  layer  is  close 
to  persistence  (no  net  change  over  the  6-h  forecast),  and  each  ozone  analysis  ratchets 
the  global  averaged  ozone  closer  to  the  global  average  of  the  observations,  reducing 
the  O-Fs  to  near  zero  (Fig.  b).  In  contrast,  for  the  experiments  with  chemistry,  the 
model  moves  the  global  averaged  ozone  in  the  upper  stratosphere  closer  to  the  global 
average  of  the  model  ozone  reference  state  r0  during  each  6-h  forecast.  Thus,  with 
chemistry,  the  assimilation  system  is  able  to  maintain  a  difference  between  the  biased 
input  data  and  the  model  forecast  (Fig.  a). 

While  the  global  mean  O-Fs  without  chemistry  converge  to  the  same  near  zero 
value  regardless  of  the  data  bias  (Fig.  7b),  the  two  ozone  analyses  are  very  different. 
Figure  8  shows  difference  plots  of  the  zonal  mean  analyses  for  15  September  2002 
for  assimilations  based  on  the  two  different  SBUV/2  data  products  in  Fig.  6.  The  input 
data  bias  between  these  two  experiments  produces  only  a  small  difference  in  the  ozone 
analyses  for  GOATS  runs  with  chemistry  (Fig.  8a),  whereas  the  ozone  analyses  without 
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chemistry  incorporate  much  of  the  bias  evident  in  the  original  data  in  Fig.  . 

To  assess  changes  at  all  the  SBUV/2  layers,  hereafter  O-Fs  over  the  time  period 
10-30  September  2002  (shaded  gray  in  Fig.  )  are  averaged  and  plotted  as  vertical 
profiles.  The  vertical  profiles  are  all  normalized  by  the  same  mid-September  SBUV/2 
values  (Table  )  to  compensate  for  the  small  absolute  O-F  values  at  upper  levels. 

Figure  9  plots  vertical  profiles  of  mean  ozone  O-Fs  averaged  over  10-30  Septem¬ 
ber  2002.  When  globally  averaged  (Figs,  a  and  b)  the  O-Fs  without  chemistry  (blue 
curves)  are  close  to  zero,  particularly  in  layers  7-10,  and  independent  of  whether  the 
input  observations  are  biased  (Fig.  a)  or  not  (Fig.  b).  The  global  bias  in  the  input  data 
can  be  seen  by  examining  the  yellow  curves  in  Figs,  a  and  b,  which  show  the  differ¬ 
ence  between  the  observations  and  the  September  zonal-mean  ozone  climatology,  r0 
(see  Fig.  3a).  Since  this  ozone  climatology  is  unchanged  for  the  experiments  shown  in 
Figs,  a  and  b,  differences  in  the  yellow  curve  between  the  two  panels  reflect  the  bias 
between  the  two  SBUV/2  data  sets.  This  data  bias  difference  is  mostly  confined  to  lay¬ 
ers  7  and  above  since  the  yellow  curves  in  layers  3-6  are  almost  identical  in  Figs.  9a 
and  9b.  In  Fig.  9a  the  bias  between  the  observation  (n16v61608)  and  the  model’s 
ozone  climatology  in  layers  7-10  (yellow  curve)  is  ~  1 0—1 5%,  while  in  Fig.  9b  the  bias 
between  the  observations  (n16v61814)  and  the  model’s  ozone  climatology  in  layers 
7-10  (yellow  curve)  is  less  than  5%.  As  expected  from  Fig.  3b  the  fast  photochemical 

time  scales  r o3=-  [ d%~rL^  |  ]  at  upper  levels  relax  the  mean  O-Fs  with  chemistry 
(red  curves)  in  Figs.  9a  and  9b  toward  the  chemistry’s  reference  state  climatology  r0 
(yellow  curves)  in  the  upper  stratosphere. 

Figures  9c  and  9d  plot  the  mean  O-Fs  averaged  over  the  equatorial  region  (between 
±15°  latitude).  This  time  the  O-Fs  for  the  experiments  without  chemistry  (blue  curves) 
are  not  constrained  to  be  close  to  zero,  as  there  can  be  a  net  advective  change  in  a 
localized  region.  Note  that  both  blue  curves  (no  chemistry  run)  in  Figs.  9c  and  9d  are 
essentially  the  same,  so  that,  as  in  the  globally  averaged  case,  the  regionally  averaged 
O-Fs  are  independent  of  the  input  data  bias  since  equatorially  averaged  advection 
patterns  are  the  same  in  both  experiments.  This  makes  sense  because  the  advection- 
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only  ozone  forecast  model  depends  only  on  ozone  gradients,  not  on  the  mean  ozone 
value.  With  the  unbiased  observations,  the  O-Fs  with  chemistry  (red  curve  in  Fig.  d) 
move  close  to  zero  in  layers  8-10  and  again  asymptote  to  the  yellow  curves  at  upper 
levels  in  both  cases. 

Figure  10  uses  the  same  presentation  as  Fig.  9  for  the  southern  and  northern  polar 
regions  (upper  and  lower  panels,  respectively).  The  southern  polar  region  is  very  dy¬ 
namically  active  during  this  time  period  (  lewman  and  Nash,  20C  ),  so  it  is  not  surpris¬ 
ing  that  even  with  the  unbiased  input  data  (Fig.  10b)  the  bias  between  the  observations 
and  the  zonal-mean  reference  climatology  (yellow  curve)  is  not  close  to  zero.  Even 
so,  the  mean  O-Fs  improve  significantly  with  the  unbiased  observations.  Note  that, 
unlike  the  global  and  equatorial  averaged  cases,  the  mean  O-Fs  with  chemistry  (red 
curves)  in  the  southern  polar  region  do  not  adjust  completely  to  the  model  chemistry. 
This  is  probably  because  the  southern  polar  dynamics  are  far  from  a  zonal  mean  at  this 
time.  Surprisingly,  even  though  regionally  (not  globally)  averaged  and  at  a  dynamically 
active  time,  the  mean  O-Fs  without  chemistry  in  the  southern  polar  region  are  very 
close  to  zero,  indicating  that  the  net  advective  ozone  forcing  in  the  model  is  small  here. 
The  northern  polar  region  shows  an  example  where  adding  chemistry,  in  addition  to 
unbiased  input  data  (Fig.  Od),  significantly  improves  the  mean  O-Fs  (red  curve)  over 
the  case  without  chemistry  (blue  curve). 

3.1 .2  Standard  deviation  of  O-F 

Figure  a  showed  some  improvement  in  the  SD  O-F  when  using  parameterized  chem¬ 
istry.  Because  dynamics  are  different  at  different  latitudes,  we  expect  the  SD  O-F  will 
depend  on  latitude,  and  therefore  this  section  focuses  on  the  dependence  of  SD  O-F 
with  and  without  parameterized  ozone  photochemistry  in  equatorial  and  southern  polar 
latitudes. 

Figure  shows  the  SD  O-F  in  the  equatorial  region  using  the  unbiased  SBUV/2 
data.  We  find  that  adding  chemistry  reduces  the  SD  O-F  in  layers  8-12,  with  much 
of  this  reduction  resulting  from  photochemistry  reducing  the  forecast  model  variance: 
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This  can  be  seen  by  expanding  the  O-F  variance  as: 

(O  -  F)'2  =  O' 2  +  F'2  -  2(0' F'),  (2) 

where  the  primes  denote  fields  after  an  area  mean  has  been  subtracted  and  the  over¬ 
bar  denotes  an  area  average.  Dashed  curves  in  Fig.  1 1  show  just  the  square  root  of 

the  forecast  model  variance,  (F'2)1^2,  for  runs  with  and  without  chemistry  (red  and  blue 
dashed  curves,  respectively).  The  observed  reduction  of  model  standard  deviation  is 
equal  to  the  reduction  of  SD  O-F  at  layers  8-10,  and  this  reduction  of  model  standard 
deviation  is  even  greater  than  the  reduction  of  SD  O-F  for  layer  1 1 . 

More  detail  is  given  in  Table  2,  where  values  for  each  term  in  Eq.  (2)  are  calculated 

separately  for  layers  8-11.  The  variance,  O'2,  of  the  SBUV/2  ozone  observations, 
remains  the  same  for  the  cases  with  and  without  chemistry.  The  variance,  F'2 ,  of  the 
6-h  model  forecast  is  calculated  by  integrating  the  forecast  model  over  the  SBUV/2 
layers  and  interpolating  the  forecast  model  to  the  SBUV/2  observation  locations,  then 
calculating  the  variance  of  these  integrated  and  interpolated  values.  The  covariance, 
(O'F1),  is  calculated  using  the  observations  and  the  interpolated  and  layer-integrated 
model  forecast  values. 

The  most  striking  effect  of  adding  chemistry  to  the  assimilation  system  is  the  huge 
reduction  of  the  forecast  model  variance.  The  covariance  (which  should  be  positive 
when  the  forecast  and  observations  are  correlated)  is  actually  negative  (anticorrelated) 
without  chemistry  for  layers  8-9.  With  chemistry,  in  layers  8  and  9,  the  (positive)  co- 
variance  becomes  the  same  order  of  magnitude  as  the  forecast  model  and  observation 
variance,  thereby  reducing  via  Eq.  (2)  O-F  variance  in  those  layers. 

Also  plotted  in  Fig.  is  the  standard  deviation  of  the  observations,  (O'2)^2  (black 
curve).  As  seen  in  Fig.  (and  the  variances  in  Table  ),  the  SD  O-F  is  larger  than 
the  SD  O  in  layers  8-12  for  the  experiment  without  chemistry.  This  means  that  a  “flat” 
forecast,  one  without  variance,  would  yield  a  better  SD  O-F  (it  would  equal  SD  O  in 
that  case)  than  using  the  no  chemistry  forecasts.  For  the  chemistry  run,  the  SD  O-F 
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is  about  the  same  as  SD  O  for  layers  6-12,  with  slightly  lower  values  in  layers  8  and  9. 
However,  the  equatorial  region  does  not  feature  such  large-amplitude  planetary  scale 
zonal  ozone  asymmetries  found  in  the  winter  extratropics,  and  therefore  these  SD  O-F 
values  are  less  than  5%  for  layers  6-1 1 . 

These  differences  in  forecast  variance  can  be  seen  in  the  forecast  ozone  fields  them¬ 
selves  (Fig.  12).  Without  chemistry  (Fig.  12a),  the  ozone  field  shows  much  more 
small-scale  structure  and  higher  variability,  with  a  difference  between  the  maximum 
and  minimum  value  over  the  field  of  3.38  ppmv;  with  chemistry  (Fig.  b),  the  ozone 
field  is  noticeably  smoother  with  a  difference  between  the  maximum  and  minimum  val¬ 
ues  over  the  field  of  1 .74  ppmv,  about  half  of  the  range  without  chemistry.  In  agreement 
with  Fig.  ,  which  shows  a  decrease  of  about  a  factor  of  3  in  time-averaged  standard 
deviation  of  the  ozone  forecast  for  layer  9  when  adding  chemistry,  the  standard  devia¬ 
tion  of  ozone  along  the  equator  at  this  one  time  in  Fig.  decreases  from  0.33  ppmv  to 
0.1 1  ppmv  when  chemistry  is  added.  Note  that  the  large  amount  of  small-scale  struc¬ 
ture  in  the  forecast  ozone  fields  without  chemistry  (Fig.  2a)  is  likely  a  consequence 
of  the  GOATS  cycle  (Fig.  ),  in  which  the  model  is  continually  adjusting  to  the  6  hourly 
re-initialization. 

Figure  plots  the  SD  O-F  and  the  SD  O  for  the  southern  polar  region:  Note  the 
change  in  scale  from  Fig.  .  This  period  corresponds  to  the  Southern  Hemisphere’s 
only  observed  major  warming  and  so,  not  surprisingly,  ozone  shows  large  variance, 
especially  in  the  lowest  layers.  At  layer  3  the  SD  O-F  is  20%  compared  to  the  SD  O  of 
near  75%.  The  model  forecast  acts  to  keep  the  SD  O-F  below  the  SD  O  for  all  layers 
for  the  chemistry  experiment,  while  the  SD  O-F  of  the  experiment  with  no  chemistry 
is  greater  than  the  SD  O  only  at  layer  12.  Like  the  equatorial  region,  the  addition  of 
chemistry  lowers  the  value  of  SD  O-F  in  the  upper  layers  (10-12)  due  to  a  reduction 
of  the  standard  deviation  of  the  forecast  model  variance  (SD  F).  In  addition,  the  model 
with  chemistry  has  somewhat  smaller  SD  O-Fs  in  layers  7  and  8  due  to  higher  covari¬ 
ances  in  those  layers  (not  shown),  suggesting  that,  in  addition  to  reducing  SD  F,  adding 
chemistry  has  improved  the  forecast  to  better  match  data  variations. 
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Looking  at  the  southern  polar  region  in  Fig.  12,  it  can  be  seen  that  without  chem¬ 
istry  (Fig.  12a),  passive  advection  forms  ozone  filaments  with  sharp  gradients,  such 
as  the  one  near  240°  E  and  75°  S.  With  chemistry  (Fig.  b),  the  same  filament  is 
much  broader,  with  much  weaker  gradients  and  a  smaller  dynamic  range.  The  slightly 
smaller,  time-averaged  SD  O-Fs  with  chemistry  (Fig.  )  show  that  photochemical 
smoothing  of  filamentary  structures  such  as  these  generally  yields  better  fits  to  the 
observations  in  the  southern  polar  region,  even  during  this  dynamically  active  time. 

3.2  Adaptive  photochemistry 

Geer  et  al.  (2006b)  used  an  ozone  assimilation  system  to  investigate  how  different 
ozone  photochemistry  schemes  affected  the  quality  of  their  ozone  analyses.  In  some 
photochemistry  schemes,  they  found  systematic  differences  between  the  scheme’s 
pre-specified  zonal-mean  reference  states  in  Eq.  (  )  and  those  corresponding  zonal 
means  of  the  GCM  forecasts  and  assimilated  observations.  These  mismatches  pro¬ 
duced  significant  errors  and  biases  in  their  final  analyzed  ozone  fields.  These  prob¬ 
lems  originate  from  assimilating  unbiased  ozone  observations  using  a  biased  ozone 
forecast  model.  In  previous  sections,  we  considered  the  opposite  problem  of  assimilat¬ 
ing  biased  observations  using  a  largely  unbiased  ozone  forecast  model. 

To  remedy  their  problems,  Geer  et  al.  (2006b)  recommended  that  attention  be  paid  in 
future  work  to  making  the  reference  states  of  the  chemistry  schemes  more  consistent 
with  those  of  the  GCM  and  assimilated  observations,  and  specifically  recommended 
that  assimilation  system  use  an  ozone  reference  state  ra  based  on  the  assimilated 
ozone  data.  This  latter  recommendation  makes  intuitive  sense,  as  it  forces  the  chem¬ 
istry  in  Eq.  (  )  to  linearize  about  reference  states  specified  by  the  latest  objectively 
analyzed  observations,  rather  than  about  internal  reference  states  that  are  completely 
independent  of  (and  possibly  biased  with  respect  to)  those  observations.  Thus,  in  this 
section  we  use  GOATS  to  test  the  Geer  et  al  (2006b)  recommendation  by  allowing 
the  ozone  and  temperature  reference  states  used  by  CHEM2D-OPP  in  Eq.  (  )  to  con¬ 
tinually  adjust  to  new  values  based  on  zonal-mean  GOATS  ozone  and  temperature 
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fields  at  recently  analyzed  times.  We  refer  to  these  GOATS  experiments  as  runs  using 
“adaptive”  (as  opposed  to  “standard”  or  “nonadaptive”)  ozone  photochemistry. 

Our  initial  implementation  of  this  concept  for  testing  in  GOATS  is  as  follows.  At 
the  start  of  each  forecast,  NOGAPS-ALPHA  takes  its  initial  state  (r,  T),  regrids  it  from 
model  levels  onto  reference  pressure  levels,  then  computes  their  zonal  means  rzm 
and  Tzm.  Since  a  single  6h  analysis  field  can  contain  high-frequency  variability,  we 
average  zonal  mean  estimates  from  a  number  of  earlier  GOATS  analysis  times  using 
the  following  recursive  time-averaging  procedure,  which  blends  the  previous  values  of 
these  reference  states  with  these  latest  zonal  means: 


(3) 


(4) 


Here  cp  is  latitude,  p  is  pressure,  A t=  6  h,  t  is  the  current  analysis  time,  and  Jav  can  be 
viewed  as  (roughly)  the  number  of  At =6  h  analysis  increments  over  which  to  apply  the 
recursive  average.  For  the  runs  here,  we  chose  Jav= 4  (JavAt=  24  h).  For  a  cold  start 
(t= 0),  the  initial  condition  we  use  at  the  first  time  step  is  r0=rzm  and  T0=Tzm . 

Only  the  state  from  the  previous  6  h  analysis,  rQ(<p,  p,t  -  At),  is  read  in  and  used  by 
this  averaging  algorithm.  While  this  greatly  simplifies  the  implementation  in  NOGAPS- 
ALPHA,  the  drawback  is  that  the  time  averaging  window  is  a  decaying  exponential, 
rather  than  a  sharp  boxcar-like  function  between  times  t-(Jav  -  1)A t  and  t.  Nonethe¬ 
less,  the  value  of  Jav  does  control  the  temporal  decay  rate  of  this  exponential,  and 
hence  the  width  of  the  time  averaging  window.  Since  NOGAPS-ALPHA  saves  (and 
thus  can  restore)  both  the  zonal  mean  and  reference  state  fields  from  all  the  previous 
runs,  much  more  sophisticated  time  averaging  algorithms  for  these  adaptive  reference 
states  could  be  devised  and  implemented.  However,  to  demonstrate  and  test  the  con¬ 
cept  in  GOATS,  we  settled  on  the  simplicity  afforded  by  Eqs.  (3)  and  (  ). 

Note  that  the  photochemical  reference  states  r0  and  T0  are  now  time  averages  of 
output  fields  from  the  GOATS.  Thus,  the  photochemistry  is  slaved  to  the  GOATS  output, 
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and  the  two  are  now  tightly  coupled.  With  r0  now  set,  Z0  is  computed  at  each  model 
time  step  by  vertically  integrating  the  r0  profile  at  each  location,  and  the  pressure- 
gridded  values  of  (r0,  T0,  Z0)  are  interpolated  at  each  model  time  step  and  location  to 
the  model’s  hybrid  o-p  vertical  levels. 

Figure  shows  the  time  series  of  the  global  averaged  SBUV/2  layer  8  mean  O-F 
results  as  in  Fig.  a  with  the  addition  of  the  corresponding  results  from  the  experiment 
using  the  adaptive  ozone  and  temperature  photochemical  reference  states  (3)  and  (4). 
This  adaptive  experiment  was  run  with  the  high-biased  SBUV/2  data,  n16v61608.  The 
change  from  the  blue  curve  to  the  green  curve  shows  the  reduction  in  mean  O-F  that 
occurs  as  the  photochemistry  adapts  to  the  observations  by  adjusting  the  reference 
states.  After  two  months  of  6-hourly  assimilation  cycles,  our  test  method  of  modifying 
the  ozone  and  temperature  climatology  has  not  yet  reduced  the  global  mean  O-F’s  to 
zero.  However,  the  mean  O-Fs  may  still  be  decreasing  slowly,  and  the  value  is  less 
than  half  that  of  the  biased  data  run  without  the  adaptive  algorithm. 

Figure  15  shows  how  the  adaptive  ozone  reference  state  (3)  (red  curve)  has  changed 
from  the  value  it  has  in  the  non-adaptive  case  (yellow  curve)  by  the  end  of  the  run 
(2  November  2002).  The  adaptive  ozone  reference  state  values  are  larger  at  all  lati¬ 
tudes  and  better  match  the  latitudinal  structure  of  the  observations  (blue  curves).  The 
adaptive  ozone  reference  state  remains  slightly  smaller  than  the  observations  in  agree¬ 
ment  with  the  non-zero  mean  O-Fs  seen  in  Fig. 

These  results  show  that  the  adaptive  ozone  photochemistry,  as  initially  implemented 
here,  produces  realistic  ozone  analyses  in  GOATS.  Adaptive  photochemistry  reduces 
the  mismatches  between  ozone  observations  and  the  photochemical  reference  states 
(Fig.  ),  which  in  turn  reduce  mean  O-Fs  (Fig.  ).  Thus,  as  suggested  by 
et  al.  (2006b),  when  the  original  photochemical  reference  states  are  biased  high  or 
low,  adaptive  photochemistry  should  improve  ozone  analyses  at  higher  altitudes  when 
assimilating  largely  unbiased  ozone  (and  temperature)  data.  However,  our  results  re¬ 
veal  a  disadvantage  of  adaptive  photochemistry  when  assimilating  biased  ozone  data. 
In  this  case,  the  adaptive  photochemistry  incorporates  this  observational  bias  into  its 
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reference  states  (Fig.  ),  and  hence  into  the  ozone  forecasts,  so  that  the  final  ozone 
analysis  fields  retain  some  of  the  observational  biases  (Fig.  ).  By  contrast,  ear¬ 
lier  GOATS  experiments  that  used  an  internal  unbiased  photochemical  reference  state 
were  able  to  correct  these  observational  biases  during  the  forecast  stage  to  yield  a 
more  unbiased  ozone  analysis  (see  Figs.  7,  f  a  and  ).  Thus,  systems  continually 
assimilating  ozone  and  temperature  data  that  are  prone  to  sudden  offsets  or  slow  drifts 
may  be  better  off  using  nonadaptive  photochemistry  to  maintain  forecast  reliability,  and 
using  mean  O-Fs  to  monitor  the  system  for  these  input  data  biases  as  in  Fig. 

On  the  other  hand,  systems  assimilating  ozone  and  temperature  data  whose  precision 
and  accuracy  are  both  high  and  stable  may  benefit  by  using  adaptive  photochemistry 
schemes.  All  these  results  reinforce  the  general  point  made  by  Geer  et  al.  (2006t)  and 
others  (e.g.,  Rood,  20(  5;  ,  )  that  biases  in  either  the  ozone  observations  or 

the  ozone  forecast  model  can  affect  the  performance  of  an  ozone  assimilation  system. 

3.3  Dependence  on  upper  boundary 

In  this  section  we  examine  the  dependence  of  the  ozone  assimilation  on  the  upper 
boundary  formulation.  We  present  results  in  which  the  model  forecast  dynamics  near 
the  stratopause  were  degraded  through  changes  in  the  initial  meteorological  fields  used 
for  the  6-h  forecast. 

As  described  in  Sect.  2,  in  initializing  the  NOGAPS-ALPHA  GCM  for  each  6-h  fore¬ 
cast  we  used  analyzed  meteorological  fields  (from  NAVDAS  and  GEOS4)  up  to  a  cer¬ 
tain  level,  which  were  then  progressively  blended  with  climatology  above.  We  ran  the 
GOATS  in  two  configurations:  one  with  that  transition  at  0.2  hPa  (the  highest  pressure 
level  available  from  the  GEOS4  analysis)  and  a  second  with  a  degraded  set  of  initial 
conditions,  in  which  that  transition  altitude  was  lowered  to  1  hPa.  These  two  exper¬ 
iments  also  differed  in  the  way  the  dynamical  fields  were  initialized:  for  the  0.2  hPa 
case  the  NOGAPS-ALPHA  non-linear  normal  mode  initialization  (NNMI)  scheme  was 
used,  while  for  the  1  hPa  case,  the  NNMI  was  deactivated.  Thus,  the  run  with  the  more 
realistic  winds  near  the  stratopause  also  benefited  from  a  more  dynamically  balanced 
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initial  state  produced  by  the  NNMI.  These  initial  conditions  are  used,  not  just  at  the 
start  of  the  experiment,  but  at  the  start  of  every  6-h  assimilation  cycle.  Both  of  these 
experiments  were  run  without  chemistry  to  focus  on  changes  to  the  dynamics. 

Figure  shows  the  vertical  structure  of  the  O-Fs  averaged  from  10-30  September 
2002.  The  mean  O-F  results  for  the  0.2  hPa  interface  experiment  (red  curve,  Fig.  b) 
are  from  the  same  experiment  as  shown  in  Fig.  9b.  Using  a  more  realistic  meteo¬ 
rological  initialization  field  from  1-0.2 hPa  improves  the  global  SD  O-Fs  at  SBUV/2 
layer  8  and  above  (Fig.  i  6a).  At  layer  1 2  the  standard  deviation  is  reduced  from  almost 
80%  to  less  than  30%.  Most  of  this  reduction  is  due  to  the  reduction  in  variance  in 
the  ozone  forecast  model,  F'2 .  While  the  most  dramatic  change  brought  about  by  the 
new  initial  conditions  is  in  the  standard  deviation,  the  mean  O-Fs  (Fig.  b)  are  also 
improved  in  levels  9,  10,  and  12.  The  mean  O-F  is  very  small  in  level  8  and  about  5% 
larger  in  level  1 1 ,  where  the  sign  of  the  mean  is  changing  with  altitude.  Thus,  although 
the  changes  in  initial  meteorological  conditions  were  at  1  hPa  and  above,  changes  are 
seen  down  to  4hPa  in  the  ozone  O-Fs. 

These  results  highlight  the  importance  of  developing  accurate  3-D  meteorological 
initial  conditions  to  altitudes  well  above  the  altitude  of  interest  for  the  ozone  assimila¬ 
tion. 


4  Discussion 

The  GOATS  combines  the  GEOS  ozone  DAS,  NOGAPS/NAVDAS  meteorological  anal¬ 
yses,  and  NOGAPS-ALPHA  forecast  model  into  a  useful  prototype  system  for  testing 
several  aspects  of  ozone  data  assimilation.  Here  we  examined  the  effects  of  adding 
an  ozone  chemistry  parameterization  (both  “standard”  and  “adaptive”),  of  observa¬ 
tional  data  biases,  and  of  improving  the  meteorological  initialization  at  upper  levels. 
Both  mean  and  SD  O-Fs  were  plotted  as  functions  of  time,  SBUV/2  layers,  and  lati¬ 
tude  regions,  illustrating  the  dependence  of  the  O-F  statistics  on  assimilation  system 
changes. 
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Our  “standard”  ozone  photochemistry  parameterization  was  the  linearized  CHEM2D- 
OPP  scheme  in  which  the  reference  states  (r0,  Ta,  la)  were  specified  using  long-term 
observational  climatologies  (see  Fig.  a).  Results  show  that  adding  this  parameterized 
ozone  chemistry  to  the  assimilation  system  can  stabilize  the  ozone  DAS  at  upper  levels 
to  drifts  or  biases  in  input  data.  The  zonal  mean  assimilated  ozone  in  GOATS  runs  with 
chemistry  (Fig.  8a)  remained  relatively  unchanged  when  fed  with  either  the  biased  or 
unbiased  set  of  SBUV/2  observations.  This  is  because,  in  the  ozone  photochemistry 
parameterization,  the  second  term  on  the  right-hand-side  (RHS)  of  Eq.  (  )  is  directly 
proportional  to  ozone,  in  contrast  to  the  ozone  advection  term  which  is  proportional 
only  to  ozone  gradients.  Thus,  the  ozone  forecast  model  with  chemistry  is  able  to 
directly  force  mean  ozone  through  the  second  term  on  the  RHS  of  Eq.  (  ),  whereas  the 
advection-only  ozone  forecast  model  can  only  produce  mean  ozone  changes  through 
eddy  ozone  fluxes  or  mean  ozone  gradients,  neither  of  which  depend  on  the  overall 
mean  ozone  magnitude. 

This  difference  between  advective  and  photochemical  tendencies  can  be  especially 
dramatic  when  examining  global  mean  quantities,  where  globally-averaged  advection 
produces  only  a  globally-averaged  vertical  eddy  flux  convergence  tendency.  Advection 
only  runs  can  often  produce  relatively  small  global  mean  ozone  tendencies  (compared 
to  the  ozone  tendencies  provided  by  the  observations)  and  therefore  the  global  mean 
ozone  produced  by  the  ozone  DAS  under  these  conditions  tends  to  draw  strongly  to 
the  mean  of  the  observations  (Fig.  b).  The  corresponding  O-Fs  tend  to  be  small  in 
this  case  (Fig.  7b). 

This  leads  to  the  general  data  assimilation  question  of  how  to  interpret  the  mean 
O-Fs  of  a  state  quantity  (such  as  global  averaged  layer  ozone)  when  that  state  quan¬ 
tity’s  model  forecasted  value  is  so  tightly  coupled  to  the  observations  (as  in  the  global 
average  of  an  advection-only  ozone  model)  that  the  forecast  is  no  longer  providing  an 
independent  estimate  of  that  quantity.  Mean  O-Fs  that  are  independent  of  data  bi¬ 
ases  (for  example,  Fig.  b)  are  a  symptom  of  a  model  forecasted  value  that  is  overly 
dependent  on  the  observations.  In  this  case,  the  model  is  no  longer  adding  significant 
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information  about  the  value  of  that  quantity  to  the  data  assimilation  process,  allowing 
that  quantity’s  value  to  asymptote  to  the  observations.  In  such  a  situation  the  mean 
O-Fs  can  be  very  small  because  the  assimilation  system  lacks  independent  model 
input  on  this  quantity,  however,  these  small  mean  O-Fs  can  be  misleading  if  this  “null” 
model  result  is  misinterpreted  as  a  well-tuned  predictive  model  result. 

Adding  physics  (such  as  parameterized  chemistry)  that  can  directly  affect  the  global 
average  of  a  state  quantity  to  a  forecast  model  yields  an  independent  estimate  of  this 
global  averaged  state  quantity  to  the  DAS.  The  global  mean  O-Fs  will  often  increase 
(slightly  for  a  well-tuned  model,  see  Fig.  5b)  over  those  generated  by  a  “null”  model, 
however,  this  increase  (if  not  too  large)  is  acceptable  if  the  model  with  physics  has  im¬ 
proved  the  SD  O-F  significantly.  In  addition  the  model  with  physics  is  now  contributing 
useful  information  to  the  DAS  about  these  global  mean  quantities. 

The  situation  is  somewhat  similar  for  the  polar  and  equatorial  mean  O-Fs  in  that 
without  photochemistry  the  mean  O-Fs  are  nearly  independent  of  the  data  bias  (com¬ 
pare  the  blue  curves  in  Figs.  cand  d;  Figs,  a  and  b;  and  Figs,  cand  d).  The 
regional  mean  O-Fs  without  photochemistry  asymptote  to  the  model’s  net  advective 
flux  convergence  for  the  region.  The  ozone  tendency  from  this  regional  net  advective 
flux  convergence  balances  the  ozone  tendency  produced  by  the  mean  O-  Fs.  Note  that 
the  regional  net  advective  flux  convergence  depends  only  on  ozone  gradients  making 
it  independent  of  mean  biases.  Thus,  without  chemistry  regional  mean  O-Fs  will  only 
be  close  to  zero  for  regions  when  the  net  advective  flux  convergence  is  near  zero.  In 
the  atmosphere  these  regional  net  advective  flux  convergences  are  often  balanced  by 
photochemistry,  and  in  an  ozone  DAS  with  photochemistry  the  regional  mean  O-Fs 
can  be  near  zero  without  requiring  the  regional  net  advective  flux  convergence  to  be 
near  zero  (for  example,  Fig.  d,  in  layers  8-1 0). 

The  mean  O-F  results  reported  here  also  give  guidance  to  the  monitoring  of  an 
ozone  assimilation  system.  Without  a  chemistry  constraint,  gradual  changes  in  obser¬ 
vation  biases  cannot  be  detected  in  the  global  mean  ozone  O-Fs  (Fig.  b)  or  regional 
mean  ozone  O-Fs  (Figs,  and  )  and  these  biases  will  be  incorporated  into  the  fi- 
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nal  ozone  analysis  (Fig.  b).  Using  “standard”  linearized  photochemistry,  ozone  at  the 
levels  where  photochemistry  is  important  can  show  the  development  of  ozone  obser¬ 
vation  biases  in  the  global  mean  O-F  (Fig.  a)  while  at  the  same  time  maintaining 
more  stable  well-performed  ozone  analysis  fields  (Fig.  a).  Experiments  with  a  pre¬ 
liminary  “adaptive”  linearized  photochemistry  fall  between  these  two  limits  (Fig.  14), 
showing  some  chemical  change  but  building  observational  bias  into  the  photochemical 
reference  state  (Fig.  )  and  thus  into  the  ozone  analysis  fields. 

Of  course,  with  or  without  chemistry  the  ozone  DAS  will  show  a  jump  in  mean  O-Fs 
if  the  input  observations  change  suddenly  (e.g.,  a  step  function  jump  with  chemistry;  a 
spike  followed  by  an  exponential  decay  without  chemistry).  This  agrees  with  the  global 
mean  O-F  time  series  shown  in  Stajner  et  al.  (200^ )  for  SBUV/2  layer  5,  where  a  spike 
followed  by  a  decay  occurred  coincident  with  a  SBUV/2  instrument  calibration  change. 
In  Stajner  et  al.  (200^ )  both  the  experiments  with  and  without  chemistry  showed  a  spike 
followed  by  a  decay,  presumably  because  the  chemistry  parameterization  used  does 
not  play  a  large  role  in  layer  5  (~32-16hPa),  where  photochemical  timescales  are  long 
(see,  e.g.,  Fig.  b).  While  the  best  way  to  monitor  input  data  biases  is  to  look  at  off¬ 
line  observation  minus  climatology  statistics  (like  the  yellow  curves  in  Figs,  and  ), 
when  the  climatological  statistic  can  be  built  into  the  assimilation  cycle  through  the 
standard  chemistry  parameterization’s  reference  state,  the  assimilation’s  operationally 
generated  O-Fs  can  provide  a  convenient  and  reliable  monitoring  proxy  at  levels  where 
photochemical  time  scales  are  short  (the  red  curves  in  Figs,  and  ). 

Our  GOATS  simulations  have  highlighted  how  the  CHEM2D-OPP  standard  photo¬ 
chemistry  parameterization  in  Eq.  (  )  yields  an  ozone  assimilation  that  is  more  robust 
to  any  sudden  anomalies,  drifts  or  biases  in  the  input  observations  (see,  e.g.,  Fig.  ). 
This  useful  property  of  the  parameterized  chemistry  in  our  GOATS  runs  was  aided  by 
using  a  reference  ozone  state  r0  (Fig.  3a)  that  was  based  on  long-term  SBUV  ozone 
measurements  (Fortuin  and  Kelder,  1998),  and  thus  provided  a  reliable  measure  of 
typical  SBUV  ozone  values  anticipated  within  each  layer  (see  Fig.  ).  Thus,  in  addition 
to  using  a  chemistry  scheme  like  CHEM2D-OPP  with  reliable  photochemical  inputs, 
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our  work  has  shown  that  the  choice  of  its  reference  state  climatologies  r0  and  T0  is  also 
important.  Ideally,  the  ozone  reference  state  r0  should  be  a  long-term  observational 
mean  of  all  the  different  types  of  ozone  observations  currently  entering  the  assimila¬ 
tion,  so  that  it  incorporates  to  some  extent  the  typical  vertical  and  horizontal  resolution 
of  the  input  ozone  observations.  For  example,  if  future  GOATS  experiments  were  to 
assimilate  higher  vertical-resolution  ozone  limb  data  such  as  from  HALOE,  MLS,  or 
MIPAS  (e.g.,  Wargan  et  al ,  2005),  the  GOATS  assimilation  of  these  data  would  likely 
work  better  by  replacing  the  SBUV-based  reference  state  used  here  with  a  new  refer¬ 
ence  climatology  r0  based  on  long-term  data  from  these  instruments  (e.g., 

Russell,  2005). 

We  have  also  used  the  GOATS  to  test  the  ability  of  “adaptive”  linearized  ozone  pho¬ 
tochemistry,  in  which  r0  and  T0  adapt  over  time  to  reflect  zonal  means  of  ozone  and 
temperature  from  GOATS  runs  at  recent  times,  to  reduce  biases  in  the  mean  ozone 
O-Fs.  Results  show  a  significant  reduction  in  the  global  mean  O-Fs  (Fig.  )  and 
a  stable  ozone  reference  state  closer  to  the  observations  (Fig.  15).  Thus,  the  use  of 
adaptive  photochemistry  climatologies  is  a  viable  option  for  improving  the  mean  as¬ 
similation  at  the  expense  of  increased  sensitivity  of  the  ozone  assimilation  to  data  drift 
coupled  with  decreased  sensitivity  of  the  mean  O-Fs  to  data  drift.  An  adaptive  formu¬ 
lation  may  be  particularly  useful  when  assimilating  new  instrument  observations  with 
high  precession  and  accuracy,  minimal  drifts,  and  for  which  a  long-term  climatology 
has  yet  to  be  formulated.  The  ozone  assimilation  with  adaptive  photochemistry  will 
then  reflect  the  instrument  biases,  if  any. 

A  possible  problem  with  the  adaptive  formulation  is  that  the  photochemistry  clima¬ 
tologies  will  depend  on  the  both  the  forecast  model  dynamics  and  the  observations. 
In  Fig.  there  is  an  increase  in  latitudinal  structure  of  the  adaptive  ozone  reference 
state  compared  with  the  non-adaptive  reference  state  and  the  observations,  that  is 
likely  a  consequence  of  the  model  dynamics.  In  the  example  in  Fig.  ,  the  additional 
latitudinal  structure  in  the  adaptive  ozone  reference  state  is  not  large  enough  to  be 
a  problem  and  the  adaptive  ozone  reference  state  at  all  other  levels  was  found  to  be 
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stable  (not  shown).  However,  at  upper  levels  especially,  the  model’s  artificial  upper 
boundary  condition  may  create  an  unstable  or  unrealistic  climatology,  especially  if  the 
data  uncertainty  is  large  in  those  regions. 

The  SD  O-  Fs  with  and  without  chemistry  show  that  chemistry  can  improve  the  ozone 
forecast  model  in  two  ways: 

1.  by  substantially  reducing  the  ozone  forecast  variance  in  models  where  the  ozone 
forecast  variance  is  likely  too  high  (Figs.  and  13) 

2.  by  increasing  the  covariance  between  the  ozone  forecast  and  the  observations 
(Table  2). 

The  high  variance  in  ozone  seen  without  chemistry  is  due  to  vertical  advection  by  re¬ 
solved  gravity  waves  (see  Fig.  12a,  especially  in  the  tropics)  and  horizontal  advection 
by  planetary-scale  waves  (seen  in  the  southern  polar  region  in  Fig.  12a).  Chemistry  in 
the  upper  stratosphere  reduces  these  advectively  created  gradients  via  photochemical 
relaxation,  producing  better  agreement  with  the  observations.  For  monitoring  and  anal¬ 
ysis  it  may  be  useful  to  plot  or  state  the  value  of  the  observation  standard  deviation  as 
a  sensible  upper  limit  to  the  desirable  SD  O-F.  If  the  SD  O-F  is  greater  than  the  SD  O, 
then  the  ozone  DAS  can  probably  be  improved  (at  least  in  the  upper  stratosphere)  with 
more  forecast  model  diffusion  to  damp  spurious  small-scale  dynamics  (see  below). 

GOATS  runs  were  also  used  to  investigate  the  role  of  upper-level  meteorological  ini¬ 
tialization  on  ozone  assimilation.  Since  GOATS  re-initializes  every  6h,  blending  the 
operational  NOGAPS/NAVDAS  and  GEOS4  upper-level  meteorological  analyses  with 
wind  and  temperature  climatologies  above,  the  potential  for  dynamical  imbalances  and 
poor  or  excessively  noisy  model  transport  exists.  Such  problems  are  ameliorated  by 
making  the  transition  to  climatology  as  high  as  possible  and  performing  nonlinear  nor¬ 
mal  mode  initialization  to  enforce  dynamical  balance  and  minimize  spurious  generation 
of  gravity  waves  which  add  “noise”  to  the  forecast  (Fig.  ).  However,  there  is  some 
evidence  that  the  dynamical  adjustment  to  the  upper  level  wind  and  temperature  cli¬ 
matologies  may  still  be  degrading  the  GOATS  results  at  the  upper  levels.  For  example, 
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the  large  SD  O-Fs  are  seen  at  upper  levels  when  GOATS  is  run  without  chemistry 
(Figs.  and  ),  as  well  as  the  somewhat  anomalous  mean  O-Fs  found  in  layer  12 
(Figs,  and  ).  Diurnal  photochemistry,  which  CHEM2D-OPP  does  not  currently 
contain,  may  also  affect  layer  12. 


5  Conclusions 

GOATS  has  provided  an  opportunity  to  study  the  basic  features  of  ozone  data  as¬ 
similation  and  the  properties  of  the  O-Fs  it  generates.  We  plan  to  use  this  experi¬ 
ence  with  GOATS  as  guidance  in  developing  ozone  data  assimilation  in  a  NOGAPS- 
ALPHA/NAVDAS  environment.  The  O-F  residuals,  while  providing  convenient  and 
very  useful  information,  need  to  be  considered  as  a  part  of  a  more  comprehensive 
evaluation  of  an  assimilation  system  and  its  components  based  on  comparisons  with 
high  quality  independent  observations.  While  GOATS  is  concerned  with  ozone  as¬ 
similation  and  photochemistry,  the  results  presented  here  apply  quite  generally  to  any 
constituent  assimilation  in  which  chemistry  or  microphysics  is  added  to  an  advection 
model. 

In  summary: 

1 .  Zonal  mean  ozone  analyses  are  more  independent  of  observation  biases  and 
drifts  when  using  an  ozone  photochemistry  parameterization  (OPP). 

2.  Mean  ozone  O-Fs  are  more  sensitive  to  observation  drifts  when  using  an  OPP. 

3.  Ozone  SD  O-Fs  are  reduced  in  the  upper  stratosphere  when  using  an  OPP  due 
to  a  reduction  of  forecast  model  noise  and  to  increased  covariance  between  the 
forecast  model  and  observations. 

4.  Changing  the  OPP  reference  state  to  match  the  observations  by  using  an  “adap¬ 
tive”  OPP  scheme  reduces  the  mean  ozone  O-Fs  at  the  expense  of  zonal  mean 
ozone  analyses  being  more  susceptible  to  data  biases  and  drifts. 
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5.  The  upper  boundary  of  the  ozone  DAS  can  affect  the  quality  of  the  ozone  analysis 
and  therefore  should  be  placed  well  above  (at  least  a  scale  height)  the  region  of 
interest. 


Appendix  A 

5 

Conversion  from  GEOS4  to  NOGAPS-ALPHA  o-p  Levels 

The  GEOS4  uses  an  L55  formulation  on  a  hybrid  o-p  coordinate  q'  of  the  form 

PK+M^  =  aK+M2^  +  bK+M2^' 

Pk+t/2  =  ak+t/2  +  bk+-\/2PS’  (A1) 

10  where  k= 0 ...  55,  pk+^/ 2  are  the  interface  (half)  level  pressures,  and  ps  is  the  surface 
(terrain)  pressure.  These  half  levels  define  the  /r=1  . . .  55  model  layers  (full  levels). 
The  NOGAPS-ALPHA  hybrid  coordinate  q  takes  the  slightly  modified  form 

P(Pk+ 1/2)  =  A(qk+ 1/2)  +  B(qk+-[/2)  [ps  -Ptop]  . 

P/r+1/2  =  Ak+ 1/2  +  Bk+-\/2  [Ps  -  Ptop]  >  (A2) 

15  where  Pjop=Pi/2 's  the  top  model  level,  here  set  to  0.1  hPa. 

In  order  to  take  the  aA+1 ,2  and  bk+y 2  coefficients  from  the  GEOS4  L55  model  and 
use  them  to  produce  the  same  interface  pressures  in  NOGAPS-ALPHA  for  a  given  ps, 
we  equate  (  )  and  (  ).  This  yields 

®/f+i/2  =  bk+ 1/2.  (A3) 

20  Ak+i/ 2  =  ak+-\/ 2  +  bk+-i/ 2  Ptop-  (A4) 

Thus,  to  produce  the  GEOS4  L55  levels  in  NOGAPS-ALPHA,  we  read  in  the  lookup 
table  of  GEOS4  ak+^/2  and  bk+-[/2  coefficients  (see  Table  3.1  of  Bloom  et  al.,  2005)  and 
convert  them  into  corresponding  Ak+y2  and  Bk+^2  coefficients  using  (A3)  and  (A4). 
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Table  1.  Mid-September  2002  values  of  SBUV/2  observations  for  layers  3-12  used  in  scaling 
vertical  profile  plots. 


Layer 

Number 

Pressure 
Range  (hPa) 

Typical 
Value  (DU) 

3 

63.3-127 

45 

4 

31 .7-63.3 

55 

5 

15.8-31.7 

65 

6 

7.92-15.8 

46 

7 

3.96-7.92 

24 

8 

1 .98-3.96 

10 

9 

0.99-1.98 

3.3 

10 

0.495-0.990 

0.96 

11 

0.247-0.495 

0.29 

12 

0.0-0.2467 

0.12 

1100 


Table  2.  10-30  September  2002  equatorial  (±15°  latitude)  variances:  ( O  -  F)'2,  F'2,  O'2,  and 
20' F1.  Values  in  DU2  are  scaled  by  104  for  layer  11,  103  for  layer  10,  and  102  for  layers  8  and  9. 


Exp 

Layer 

(O  -  F)'z 

20' F' 

F'2 

O'2 

No  Chem 

11 

17.79 

1.81 

18.70 

0.90 

Chem 

11 

1.14 

-0.05 

0.19 

0.90 

No  Chem 

10 

18.45 

-0.51 

17.74 

0.20 

Chem 

10 

0.36 
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Fig.  1.  Schematic  depiction  of  a  GOATS  cycle.  The  two  main  components  are  NOGAPS- 
ALPHA  (ozone  forecast)  and  GEOS  ozone  DAS  (ozone  analysis).  SBUV/2  layer  ozone  obser¬ 
vations  are  assimilated  into  the  cycle.  The  ozone  fields  continue  through  the  cycle,  however, 
the  meteorological  fields  are  re-initialized  every  6  h. 
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(a)  NOGAPS:  L30  a  levels  (b)  NOGAPS-ALPHA:  L55  hybrid  o-p  levels  (c)  Pressure  Height  Layer  Thicknesses  (d)  December  Ozone  Climatology  (ppmv) 
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Fig.  2.  Model  half  levels  around  34.5°  N  latitude  circle  for  (a)  NOGAPS  L30  (pure  a  levels) 
and  (b)  NOGAPS-ALPHA  L55  levels  (hybrid  o-p  levels).  Orange  levels  demark  heavily  dif¬ 
fused  upper  layers,  red  line  shows  top  half-level  pressure  pT0P,  yellow  line  in  (b)  shows  lowest 
isobaric  half  level  (177hPa).  (c)  Model  layer  pressure-height  thicknesses  A Zk  for  these  model 
formulations,  as  well  as  L60  NOGAPS-ALPHA  levels,  (d)  Zonal-mean  daytime  ozone  volume 
mixing  ratio  climatology  (ppmv)  for  December  as  function  of  latitude  and  pressure. 
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(b)  CHEM2D-OPP  t03  (days):  September 


Fig.  3.  (a)  Ozone  reference  state,  r0  (ppmv)  and  (b)  ozone  photochemical  relaxation  time,  r0 
(days)  as  function  of  latitude  and  pressure  for  September.  The  gray  bands  denote  the  SBUV/2 
layers  (see  Table  for  details). 


1104 


^■n  i  i  i  1 1  n  :  i  i 

150  190  230  270  310  350  390  430 

TOZ  [DU] 

Fig.  4.  Total  ozone  in  the  Southern  Hemisphere  for  20-30  September  2002  (48-h  intervals) 
from  the  Total  Ozone  Mapping  Spectrometer  (TOMS;  top  row)  and  from  the  GOATS  T79L55 
run  with  CHEM2D-OPP  chemistry  (bottom  row). 
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Fig.  5.  Comparisons  of  SBUV/2  layer  8,  global  averaged  (a)  standard  deviation  and  (b)  mean 
O-Fs  (DU)  for  GOATS  without  (blue  curve)  and  with  (red  curve)  an  ozone  photochemistry  pa¬ 
rameterization.  A  2-day  boxcar  filter  has  been  applied  to  the  6-hourly  global  average  O-Fs  to 
reduce  diurnal  fluctuations,  and  the  initial  day  is  not  plotted.  The  gray  shaded  regions  denote 
the  domain  that  is  time  averaged  in  later  figures.  SBUV/2  layer  8  represents  the  vertically  inte¬ 
grated  ozone  from  3.96-1 .98  hPa,  with  an  typical  global  average  September  value  of  ~10  DU. 
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Fig.  6.  SBUV/2  layer  8  ozone  observations  (DU)  on  15  September  2002  plotted  as  a  function 
of  latitude  for  versions  n16v61608  (blue  curves)  and  version  n16v61814  (red  curves).  Also 
plotted  are  corresponding,  vertically-integrated  over  SBUV/2  layer  8,  September  background 
ozone  values,  r0,  based  on  the  (  )  zonal-mean  ozone  climatology  (yellow 

curve)  used  in  the  CHEM2D-OPP  scheme  in  Eg.  (  );  see  Fig.  a. 
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Fig.  7.  Comparisons  of  SBUV/2  layer  8,  global  averaged  mean  O-Fs  (DU)  for  GOATS  runs 
(a)  with  and  (b)  without  the  ozone  photochemistry  parameterization  (CHEM2D-OPP).  The  blue 
curves  assimilated  the  high-biased  (n16v61608)  version  of  SBUV/2  observations.  The  red 
curves  assimilated  the  more  recent  SBUV/2  n16v61814  observations.  A  2-day  boxcar  filter 
has  been  applied  to  the  6-hourly  global  average  O-Fs  to  reduce  diurnal  fluctuations,  and  the 
initial  day  is  not  plotted.  The  gray  shaded  regions  denote  the  domain  that  is  time  averaged 
and  plotted  in  other  figures.  SBUV/2  layer  8  represents  the  vertically  integrated  ozone  from 
3.96-1 ,98hPa,  with  a  typical  global  average  September  value  of  ~10  DU. 
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Fig.  8.  Differences  in  zonal-mean  ozone  analysis  fields  on  15  September  2002  between 
GOATS  assimilations  of  high-biased  SBUV/2  observations  (n16v61608)  and  unbiased  SBUV/2 
observations  (n16v61814)  for  runs  (a)  with  and  (b)  without  the  ozone  photochemistry  parame¬ 
terization  (CHEM2D-OPP).  The  contour  interval  is  0.1  ppmv. 
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Fig.  9.  Vertical  profiles  of  ozone  mean  O-Fs  expressed  as  a  percentage  relative  to  the  nominal 
mean  values  in  Table  ,  averaged  over  10-30  September  2002  with  (red  curves)  and  without 
(blue  curves)  parameterized  ozone  photochemistry.  Plots  (a)  and  (c)  show  results  from  as¬ 
similating  the  earlier  (high  biased)  version  of  SBUV/2  observations  (n16v61608);  plots  (b)  and 
(d)  show  results  from  assimilating  the  more  recent  (unbiased)  version  of  SBUV/2  (n16v61814). 
Plots  (a)  and  (b)  are  global  averages,  plots  (c)  and  (d)  are  equatorial  averages  (15°S-15°N). 
Also  plotted  is  the  mean  difference  between  SBUV/2  observations  and  the 
(  )  ozone  climatology  (yellow  curve)  for  the  corresponding  data  sets  and  regions.  The 

gray  horizontal  lines  mark  the  boundaries  of  the  SBUV/2  layers  (Table  )  labeled  by  the  gray 
numbers  on  the  right  of  each  plot.  inn 


Fig.  10.  Vertical  profiles  of  ozone  mean  O-Fs  averaged  over  10-30  September  2002  with  (red 
curves)  and  without  (blue  curves)  parameterized  ozone  photochemistry.  Plots  (a)  and  (c)  show 
results  from  assimilating  the  earlier  (high  biased)  version  of  SBUV/2  observations  (n16v61608); 
plots  (b)  and  (d)  show  results  from  assimilating  the  more  recent  (unbiased)  version  of  SBUV/2 
(n16v61814).  Plots  (a)  and  (b)  are  averaged  from  90°S-60°S,  plots  (c)  and  (d)  are  averaged 
from  60°  N-90°  N.  Also  plotted  is  the  mean  difference  between  SBUV/2  observations  and  the 
(  )  zonal-mean  ozone  climatology  (yellow  curve)  for  the  corresponding 

data  sets  and  regions.  The  gray  horizontal  lines  mark  the  boundaries  of  the  SBUV/2  layers 
(Table  )  labeled  by  the  gray  numbers  on  the  right  of  each  plot. 
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Fig.  11.  Vertical  profiles  of  the  standard  deviation  of  ozone  O-Fs,  [( O  -  F)'2]1/2,  averaged  over 
10-30  September  2002  and  within  ±15°  latitude  with  (red  curve)  and  without  (blue  curve)  pa¬ 
rameterized  ozone  photochemistry.  Both  curves  result  from  assimilating  the  unbiased  SBUV/2 
observations.  The  corresponding  colored  dashed  lines  plot  the  standard  deviation  of  the  ozone 

forecast  interpolated  to  the  observation  locations,  (F'2)1^2.  The  dashed  lines  are  often  close  to 
the  corresponding  solid  [( O  -  F)'2]1^2  curve.  The  black  curve  plots  the  standard  deviation  of  the 
SBUV/2  observations,  (O'2)1/2.  The  solid  red  curve  lies  beneath  the  black  curve  at  layers  11 
and  12.  The  gray  horizontal  lines  mark  the  boundaries  of  the  SBUV/2  layers  (Table  )  labeled 
by  the  gray  numbers  on  the  right  of  each  plot. 
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Fig.  12.  Ozone  (ppmv)  at  1 .65  hPa  taken  from  a  6-h  forecast  valid  on  24  September  2002  for 
GOATS  experiments  run  (a)  without  and  (b)  with  an  ozone  photochemistry  parameterization. 
The  contour  interval  is  0.25  pppmv.  Darker  contour  curves  are  drawn  at  4,  5,  and  6  ppmv. 
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Fig.  13.  Vertical  profiles  of  the  standard  deviation  of  ozone  O-Fs  averaged  over  IQ- 
30  September  2002  and  from  90°  S-60°  S  with  (red  curve)  and  without  (blue  curve)  param¬ 
eterized  ozone  photochemistry.  Both  curves  are  based  on  the  unbiased  SBUV/2  observations. 
The  black  curve  plots  the  standard  deviation  of  the  SBUV/2  observations,  [0,z]1/2.  Note  the 
change  in  x-axis  scale  from  Fig.  .  All  curves  are  normalized  by  the  layer  values  in  Table  . 
The  gray  horizontal  lines  mark  the  boundaries  of  the  SBUV/2  layers  (Table  )  labeled  by  the 
gray  numbers  on  the  right  of  each  plot. 
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Fig.  14.  Same  as  Fig.  a  with  the  addition  of  results  taken  from  a  GOATS  run  using  CHEM2D- 
OPP  modified  to  incorporate  adaptive  photochemical  reference  states  for  ozone  and  tempera¬ 
ture  (green  curve). 


1115 


Latitude 


Fig.  15.  SBUV/2  layer  8  ozone  observations  (DU)  on  2  November  2002  plotted  as  a  function 
of  latitude  for  version  n16v61608  (blue  curves).  Also  plotted  are  the  corresponding,  vertically- 
integrated  over  SBUV/2  layer  8,  background  ozone  values,  ra,  based  on  the 
(  ),  zonal-mean  ozone  climatology  (yellow  curve),  and  the  2  November  2002  adaptive  ref¬ 

erence  state  (red  curve)  used  to  specify  the  CHEM2D-OPP  relaxation  term  in  Eg.  (  ). 
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Fig.  16.  Global  (a)  standard  deviation  and  (b)  mean  ozone  O-Fs  averaged  from  IQ- 
30  September  2002  normalized  by  the  average  SBUV/2  ozone  for  mid-September  2002  (Ta¬ 
ble  )  and  expressed  as  a  percent,  plotted  as  a  function  of  pressure.  Both  experiments  are 
based  on  the  unbiased  SBUV/2  observations  (n16v61814).  The  gray  horizontal  lines  mark  the 
boundaries  of  the  SBUV/2  layers  (Table  )  labeled  by  the  gray  numbers  on  the  right  of  each 
plot. 
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